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The interpretation of optical conductivity in the normal and superconducting states 
is considered in the frame of the standard Isotropic Single wide Band (ISB) model using 
the theory proposed by Nam (Nam S.B., Phys. Rev. 156 470-486 (1967)). The exact 
analytical inversion of the normal state Nam equations is performed and applied for the 
recovery of the reliable information from the experimental data. The Allen formula is 
derived in the strong coupling approximation and used for the physically transparent in- 
terpretation of the FIR absorption. The phenomenological "generalised" Drude formula 
is obtained from the Nam theory in a high temperature approximation. It is shown, that 
the reconstruction of the shape of the spectral function a 2 F(u>) from the normal state 
optical data at T>T C is not unique and the same data can be fitted by many spectral 
functions. This problem is considered in detail from different points of view. At the same 
time, using the exact analytical solution, one can get from the normal state data a useful 
piece of information, namely, the value of the coupling constant, the upper energy bound 
of the electron-boson interaction function, and the averaged boson frequency. Moreover, 
they are even visually accessible, if the optical data are presented in term of the frequency 
dependent optical mass and scattering rate. The superconducting state Nam formalism 
and related simplified theory are analysed from the user point of view. A novel adaptive 
method of the exact numerical inversion of the superconducting state Nam equations 
is presented. Since this approach uses the first derivative from the experimental curve, 
the signal/noise ratio problem is discussed in detail. It is shown that, the fine structure 
of the spectral function can be recovered from optical data in the case of s-wave pair- 
ing symmetry. In contrast, in the d-wave case the resulting image is approximately the 
convolution of the input electron-boson interaction function and the Gauss distribution 
exp(— uj' 2 /Ao pt ). A simplified visual accessibility (VA) procedure is proposed for "by eye" 
analyses of the superconducting state optical reflectivity of the ISB metals with s and 
d pairing symmetry. The bosons, responsible for the superconductivity in YBa2Cua07, 
exhibit a "phonon-like" spectral function with the upper frequency bound less than 500 
cm -1 and the averaged boson frequency near 300 cm -1 . 



I. INTRODUCTORY REMARKS, RESERVATIONS AND BASIC EQUATIONS 

• This is the preliminary version of invited lecture for Nato-ASI "Material Science, Fundamental 
Properties and Future Electronic Applications of High-Tc Superconductors", Albena, Bulgaria, 
September 1998, Kluwer Academic Publishers, Dortrecht, pp. 323-360 (2001) 

• The interpretation of the normal and superconducting state Far Infrared Region (FIR) properties 
of metals is still rather ambiguous. In this lecture I will restrict myself to the consideration of the 
of Isotropic Single wide Band (ISB) model based on the Nam formalism [1]. 

• From the optical conductivity o{uj) one derives the so called transport spectral function a 2 r F{u>). 
It has the same spectral structure, i.e. the same number of peaks with similar relative positions 
(see Fig. 4 in [2] and Fig.l in [3]) as a standard a 2 F{u>) which enters the Eliashberg equations. But 
their amplitudes can be lower, and as a result the transport coupling constant Aj r is less than the 
standard A. In this lecture I will not distinct the a 2 F(oj) and a 2 r F(ui). 

• If the electronic band is very narrow, the forward scattering could dominate. In this case the optical 
properties manifests themselves by the negligible coupling, while the tunnelling spectroscopy gives 
reasonable finite coupling strength. A similar problem arises when the quasiparticles from the flat 
or nested parts of the Fermi surface are mainly responsible for the superconductivity. Since they 
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have small Fermi velocities, their influence in the transport response function which is proportional 
to v F is weak. Such a scenario was recently proposed for borocarbides [4] . 

• In the frame of the Eliashberg theory it is impossible to elucidate the nature of the bosons which 
are responsible for the superconductivity. In this context I will use the term boson, keeping for 
short the standard electron-phonon notation a 2 F(u). 

• I will restrict myself to the consideration of the local (London) limit case, that is, to the normal skin 
effect scenario. The generalisation of results obtained to the nonlocal (Pippard) limit is possible 
but will not be considered here due to the lack of space. 

• The power of the proposed methods is illustrated analysing available experimental data for a single- 
domain YBa2Cu307_(j crystal [5]. Note, that to the best of my knowledge, the experimental HTSC 
superconducting state reflectivity was not been considered from the inverse problem viewpoint at 
all, despite the single early attempt [6] . 

• Theoretical analysis of any physical quantity, as optical conductivity, is performed in frame of some 
model which includes a theoretical formalism and some assumptions about the band structure, the 
anisotropy of the coupling function and so on. The considered here ISB model is the standard 
default. Nevertheless, in frame of its s-wave version the question: "What is the value of the gap? 
" has an unique answer, but for <i-wave case the correct questions are either :"What are the values 
of the gaps? " or "What is the value of the maximum gap? ". In other words, the absolute values 
and the physical meaning of the quantities are model dependent. 

• If the model is reach enough, one can calculate many physical properties using the same set of 
input (material) parameters. A solution of the inverse problem means the determination of these 
material parameters from the experimental data. A short list of the ISB model input parameters 
is given in section I A. 

• As a rule, the inverse problem is ill-posed. The experimental data contain less information, than 
we would like to know. The reasons for are the temperature induced broadening and the finite 
signal to noise ratio S/N. 

• In this lecture we will mainly be interested in the amount of the remaining information about 
the spectral function a 2 F(uj) 1 which can be obtained from the optical data. In spectroscopy the 
unit of information is a peak. A single peak needs four parameters for its description: the upper 
and the low frequency bounds fi mnl and 0, m in, the position of the maximum and its amplitude. 
An narrow symmetric peak within the broad band needs three parameters: the amplitude, the 
halfwidth, and the position. The sum of all peaks yields the upper and lower frequency bounds 
Qmax and ^l m in of the spectral band. Since the fl m in of the a 2 F(uj) is usually equal to zero, the 
remaining three band parameters of a 2 F(ui) can be determined even from the normal state data 
[7]- 

A. Isotropic single band model 

The standard ISB model [8] is the most developed part of the modern theory of superconductivity. It 
describes quantitatively the renormalization of the physical properties of metals due to the electron-boson 
interaction. The input quantities of the ISB model are the density of states at the Fermi energy N(0), the 
Fermi velocity vf, the impurity scattering rate jimp, the Coulomb pscudopotential fi*, and the spectral 
function a 2 F(uj) of the electron-boson interaction. The scale of transport and optical properties is fixed 
by the plasma frequency Uuj p i = y / 4'Ke 2 v' F N(0)/3. In this section I present only the final expressions of 
the Nam theory, more detailed information can be found in [1,2,9]. 

In order to calculate any physical property in the superconducting state one needs a solutions of 
Eliashberg equations (EE) [11]. In my opinion, the following EE representation is most suitable for 
numerical solution by iterations 
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where A(w) and £>(w) are the renormalized gap function and the renormalized frequency respectively, and 
7imp denotes the impurity scattering rate within the Born approximation. The real and the imaginary 
parts of the Eliashberg functions A(w) and Q(u>) are connected by the Kramers-Kronig relations. Hence, 
they have the same Fourier images. This reasoning yields the fast solution procedure. The convolution 
type integrals (1-2) should be calculated by the Fast Fourier Transform (FFT) algorithm. The inverse 
complex Fourier transformations of the results obtained give complex values of A(ui) and Q(uj). If one 
would like to use this efficient FFT method also for d-wave pairing symmetry, the process of arithmetic- 
geometric mean [12] for the evaluation of the complex elliptic integrals ( see below Eqs. (67,68) is 
strongly recommended. The density of states ReN(oj) and the density of pairs ReD(u>) 
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could be approximated by step functions. The convolution of any spectrum with the step function 
sign(y 2 — Aq) results in the shift by A at y > and by — A at y < and its integration over 
y. It is important, that due to the presence of singularities in N(ui) and D(w) at uj — Ao the more 
complicated functions such as / duj' N(oj')N(uj — oj') and / du>' D{oj')D(uj — u>') still contain a jump (now 
at u> = 2Ao). Since the convolution is a quite general property of the Green functions approach, one can 
expect that the first derivatives from a spectrum (of different nature) over frequency will reproduce the 
input spectral function a 2 F(ui) (may be with the appropriate phase distortion). When applied to the 
optical conductivity, this approach gives visual accessibility (VA) procedure discussed in the sec. III. 
In the normal state the set (2) is reduced to the following simple formula 
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which however looks unwieldy in comparison with the same formula presented in Matsubara formalism 

u)(iu> n ) = u n = Lu n -i + 2irT(l + X n ) (6) 

where 



\ k = a 2 F{iv k ) = 2 j dza 2 F(z)z/(z 2 + v 2 k ) 



(7) 



are kernels of the spectral function, uok = 7rT(2fc + 1) and Uk = 2irTk are fermion M) and boson (i/^) 
Matsubara energies, respectively, and u>o = Jimp/2 + ttT(1 + Ao). 
The normal state ISB optical conductivity takes the following forms 
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and 

UJ 2 ™ _1 | 

cr(iu n ,T) = a n =-r E - V - — — , (9) 



n — k—1 



in the real and imaginary axes techniques, respectively. 

The tedious expression for the superconducting state ISB optical conductivity has the form [1,9,10] 
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are coherent factors. 
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g rr = l- N(x)N(x + u) - D(x)D(x + uj), 

g ra = 1 + N*(x)N(x + lo) + D*(x)D(x + u) (11) 



II. NORMAL STATE OPTICAL CONDUCTIVITY 

A. The exact solution of an inverse problem 

It is amusing, that the Eq. (9) is solvable for ujk and Eq. (6) for Afe. Hence, at least theoretically, an 
exact analytical solution exists. It takes three or four steps as shown below. 

I) One has to perform the analytical continuation of the conductivity a(u>) from the real energy 
axis, where we are living, to the poles of the Bose distribution function (Matsubara boson energies) 
ii/ k = 2'KiTk. This procedure is similar to the Kramers-Kronig analysis, which is nothing other than 
the analytical continuation from the real frequency axis to itself. For example, if one fits the data by a 
sum of Drude and Lorentz terms or by formula (33), one has simply to substitute the complex values 
ivk for the frequency. Note, that the non-metallic (IR active direct phonon contribution and intcrband 
transition) part of the dielectric permeability e p h(to) have to be subtracted from the total one e(uj) before 
starting of the analysis. In the genuine far infrared region one could use the real constant f M = e p h{0), 
instead of e p h(u)), but if the spectral range is wide, the subtraction of e p h(uj) is not trivial. 

II) The renormalized frequencies should be calculated from (9) as follows 

< 12 » 

u>i = — ^--cJ , (13) 
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uo 2 uj 2 n ~ 2 1 

^n-i — 7r~7 u7 , where A n = a„ --r^~ V! - — — ~ • ( 14 ) 

2ixA n Aim fr[u k + Un-k-i 

III) Inverting Eq. (6), one could evaluate the values of the spectral function A„ = a 2 F(v n ) we looked 
for 
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IV) Using, for example, Pade polynoms [14] one has to continue the electron-boson interaction function 
back to the real axis. 
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FIG. 1. Panel (a) shows the frequency dependence of the normal state reflectivity (solid line, T=100 K, 
Schutzmann J. et al., Phys. Rev., B46, 512 (1992)) in comparison with the model ones. The curves, shown by 
the dotted and the dotted-dashed lines, was calculated using "model 2" and "single peak" spectra (presented 
on the panel (b)). The dashed line is the least square fit of YBa2Cu30n.9 a-axis reflectivity by the favourite 
polynom Eq. 32 (with the fixed = 12). 



In practice, if the signal to noise ratio is not extremely large, the last step makes no sense due to the 
dramatic increase of the uncertainty in the A„ with elevation of n. The point is that the small quantities 
A„ Eq. (16) and A n Eq. (14) are the difference of large quantities. Note, that the uncertainty in the 
A„ values has to be small, not in comparison with A„ itself, but with the difference between A„ and its 
high energy asymptotic value 
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If the average boson energy has the order of 2irT, the last condition (SX n <C |A„ — A„|) may be violated 
even for n=3, since A„ and A„ have similar values. For example, at T =100 K for the spectral function, 
shown by the dotted line in Fig. 1, A„/A n =1.37, 1.12, 1.06, 1.03 for n=l, 2, 3, 4 respectively. Three 
or four A„ values only form a too small parameter set to recover correctly the shape of the spectral 
function. Actually, the possibility to obtain at least small amount of reliable information, is the additional 
advantage provided by the Matsubara solution in comparison with the real axis one where all values are 
unreliable [15]. Fig. 1 shows the good agreement between the experimental normal state 1-2-3 single 
crystal reflectivity data [5] and the calculated ones using for various spectral functions. Unfortunately, 
the lack of a unique solution promoted unreasonable speculations on the mechanisms in several novel 
superconducting systems in connection with proposed hand-made electron-boson interaction functions. 
In my opinion, twelve years after the discovery of the superconductivity in cuprates it is timely to return 
to a self-consistent description of the superconductivity in this compounds. 
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The first Matsubara value of the conductivity a(iv\){= a\) gives the value of the coupling constant 
A(= Aq), if the plasma frequency ui v i and the impurity scattering rate 7i mp are known a priori 



1_ = 8tHT(1 + Ao) + ^ limp 
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The reflectivity R(w) can also be continued to the imaginary frequency axis and included into the solution 
of the inverse problem. Following this way, one can derive from (19) the following formula 

X ° = ~ 2^T 1 + (2^T)2(coth 2 (V 1 /4)- eoo )' (20) 

where is the low-frequency value of the non-metallic (phonon and interband transition) part of 
dielectric permeability, and K„, is the electromagnetic kernel 
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The first Ai and the second A2 allow ( for example ) to estimate the average boson energy and very 
approximately the width of the spectral function. If there is the possibility to combine this incomplete 
information about a 2 F(iv) with the results of other spectral measurements, for example, tunnelling or 
neutron data, one can qualify the transport spectral function. The so called "model 2" electron-boson 
interaction function [7] was recovered assuming, that the solution has to be resemble the phonon density 
of states, measured by inelasitc neutron spectroscopy. 



B. Optical mass and impurity scattering rate 

In rcf. [7,15] we shown, that it is useful to describe the optical conductivity in terms of the so called 
"extended" Drudc formula 



/4tt ^p/4?r 
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Here the optical relaxation W(co), the optical mass m* pt (w, T)/rrib and the optical scattering rate 
7 opt (cj, T) as the complex, real and imaginary parts of an inverse normalised conductivity (J 2 i/4tto-(u>, T) 
have been introduced. This presentation is old and popular in spectroscopy, but the practical value of 
the Eq. (22) has not been exploited much. The optical relaxation will be compared with the effective 
mass mlff(u),T)/rrib and the effective scattering rate 7 e //(w,T), which by definition are the real and 
imaginary parts of the renormalized frequency 

^,T) = ^ ( "' r) + ^<"' r ) (23) 
m h 2 

When the frequency dependencies of j opt (u),T) and m*(uj 1 T)/mb are plotted, the coupling constant, 
the average boson energy and the upper boundary of spectral function are visually accessible (see Fig. 
2). The appropriate visual accessibility rules can be derived from the properties of the renormalized 
frequency ui using the Allen formula Eq. (26) [2]. Since such useful expression was originally obtained in 
the weak coupling approximation (A <C 1) only, below I rederive it for the general case, including strong 
coupling. 

For any k the values of Afe are positive and are of the same order of magnitude, except may be 
A . It means, that in the zero order the series uu n defined by Eq. (6) is an arithmetical progression. 
Consequently, according to the Gauss rule for the arithmetical series the denominators u>k + u> n -k-i in 
(9) at a given n do not depend on k, that is, they coincide. In view of Eq. (22), Eq. (9) takes the form 

W n +1 ~ LO k + U) n - k . (24) 

Hence we can expand the denominator in (9) 
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to the first order in power of + uj n -k-i — W n )/W n . After some simple algebra we arrive at the 
sought-for Allen formula 
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FIG. 2. The frequency dependencies of the normal state optical and effective masses (a) and scattering rates (b) 
calculated using the "model 2" spectral function with \—2 (shown in panel (a) for comparison). The quasiparticle 
effective mass m*j- /(w) /nib and scattering rate 7e//(w) at T =1 K (solid lines, panels (a) and (b)) reproduce 
the fine structure of the spectral function, but at T =100 K (dashed lines) this function are smeared out due 
to averaging over Fermi distribution. At lo > Q max the quasiparticles become "undressed", m*ff{u))/mb ~ 1 
and 7e//(w) « const. The optical mass mo P t(u)/mb and the scattering rate 'y opt (u>) even at T —1 K (dotted 
lines) do not reproduce the structure of the spectral function due to the energy averaging, all the more at high 
temperature ones (T =100 K, dashed-dotted lines). 



The real axis versions of Eq. (26) at T =0 look similar 
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and justify the following simple description of the interaction of light with the conducting subsystem 
of the normal state metal. When a photon with an energy w has been absorbed, two excited virtual 
quasiparticles, an "electron" and a "hole" are created. If the first particle's frequency is u>', the second 
ones frequency should be u> — lo' . The exited "electron" and "hole" relax to the Fermi level according 
to the quasiparticle laws (see Eqs. 5,6). Since lo' can varies from to lo, the optical relaxation W(lo) is 
the frequency-averaged (^ J^), double (electron+hole — > factor 2) renormalized frequency lo. 
At T = the ISB effective scattering rate j e ff(co) = 2ImLo{uj) 



f-LJ 

7e//(w) = Jimp + 2?r / dya 2 F(y) 
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(30) 
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has a clear physical meaning. The exited quasiparticles with the energy uj can relax by the emission of 
virtual bosons with the energy varying the range from to w or by the impurity scattering. 
At finite temperature the Allen formula 

If tanh(^) +tanh(^) 
W(to) = - / dz 2T ' n K -^-2u)(z) (31) 

U! J 2 

and the effective scattering rate (5) have exactly the same interpretation, since the only difference 
between them and Eqs. (27,27,29,30) is the presence of the Fermi distributions instead of step functions 
at T=0. 

Figure 2 shows the frequency dependencies of the normal state optical and effective masses and 
scattering rates calculated using the model spectral function [7] with A=2. The results are plotted for 
the low temperature T=l K and for T=100 K which is close to T c and has the order of flboson/^nT , 
where f2t, osorl =330 cm -1 is the averaged boson frequency Both m* e ^{uj) /nib and 7 e //(u>) show the 
characteristic features corresponding the peculiarities of the phonon spectrum. When the frequency ui 
exceeds the upper energy bound of the boson spectrum f2 max , the quasiparticle becomes "undressed". 
Its mass m* e f f (uo)/m b 1 and effective scattering rate -f e ff(uj) ~ const, but of course, the absolute 
value of j e ff (to — ► oo) could be huge. 

The optical mass m* opt {uj) / 'nib and scattering rate 7 op t(w) have no features at uo m tt max due to the 
energy averaging (27,28,29,31). One can only calculate the approximate value of £l max as the frequency 
where the m* opt {ui)/nib~ 1 and j opt (tu) reach the halves of their maximum magnitudes. Fig. 2 shows, that 
this my rule can not be applied even to the calculated frequency dependence of the optical scattering 
rate 7 opt (w), the more to the experimental curve, where the high energy behaviour is distorted by the 
interband transition contribution to the dielectric permeability. In contrast the frequency dependence 
of the optical mass m* pt (u>) / 'nib is suitable for the application of the simple criterion 

m*(w = n max )/m b - 1 « 0.5(m*(w = 0)/m b - 1). (32) 

Moreover its low energy value m* pt (w=0)/mjsjl + X tr itself contains a useful piece of information. Fig. 
2a shows, that the approximate values of the transport coupling constant Xt r ~ ( m opt (0) l m b — 1) and the 
upper frequency bound of the electron-boson coupling function [m* pt (Q max )/mb — 1 ~ 0.5(m* pt (0)/mb — 
1)] are visually accessible. If the spectral function is wide and its peak position (or averaged frequency) 
coincide roughly with the Q max /2, the optical relaxation W(u>) is smooth and the elucidation of Qboson at 
high temperature is impossible. When an exotic superconductor with the narrow, (5-function like spectral 
function will be discovered, following this way one can easily construct the recipe for the evaluation of 
the strength and the position of the single peak. 



C. My favourite polynom 

The Pade polynom analytical continuation [14] is a capricious and not completely correct procedure. 
Nevertheless, if the temperature is high enough in comparison with the energy under consideration it 
works splendidly. The numerical experiment has given to my great surprise the following (and to some 
extent obvious) result (see Fig. 3). The good agreement between the real axis calculations and the 
Matsubara ones takes place every time when the approximation using 1000 a n points gives the same 
result as the four points Pade approximation. 

This remarkable Pade polynom with only four parameters a^ ; /47r, A, B and C 

*M = A _t + *-^ (33) 

fits well the majority of the experimental data at T > T c and itself is the demonstration of the result 
obtained in section II A: the few first Matsubara values are responsible for the frequency dependence 
of a(u). 

The Pade polynom (33) is a formula of merit. At first, it has the correct analytical properties. Its 
two poles are located in the lower half-plane of the complex energy. The high (u> — > oo ) and low 
(w — > 0) energy asymptotic behaviour arc reasonable [23]. At second, it is ideal for the least square 
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fits of experimental R(w) [24], as well as for simple analytical estimations. Having the exact solution 
(see sec. II A) we have to substitute the complex values ivk for energy, perform easily the analytical 
continuation to the imaginary axis and analyse the values of tTj . In the way we obtain the interpretation 
of A, B, C. Another possibility is a fast and easy analysis of the data in terms of the optical mass 
m* pt (u),T)/mb and the optical scattering rate j opt (u),T) 

rKpti^T) . BC 

lopt {u,T) = A+-^^. (35) 
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FIG. 3. Panel (a) shows the frequency dependencies of the real and imaginary parts of the calculated normal 
state optical conductivity using the real axis formalism (solid lines) and Matsubara technique (squares - Recr(uj), 
circles - Imcx(w)). The analytical continuation from the imaginary axis to the real one was made by the Pade 
polynom with the degree of polynomial n=4. Calculations was made for T=100 K using "model 2" spectral 
function. Panel (b) shows the frequency dependence of the normal state reflectivity R(w) of RbsCeo measured 
at T=40 K by Degiorgi et al., Phys. Rev. Lett. 69, 2987 (1992) (circles ) in comparison with the least square fit 
by the favourite polynom. Since four parameters were too much, the value of the plasma frequency uj v i =10000 
cm -1 was fixed. The obtained frequency dependence of the optical mass is shown in the inset. 



At third, if the low-frequency reflectivity R(w) measured at high temperature can not be fitted by Eq. 
(33) or the obtained value of Ao contradicts the generally accepted values, it might indicate, that the 
material under consideration is not a standard Isotropic Single wide Band (ISB) metal and its properties 
have to be analysed if terms of another, unconventional model. In order to illustrate the utility of the this 
approach, let us consider another superconducting compound with relatively high transition temperature 
Ro?,Cqq [16,18,17]. Since the ratio of the experimental reflectivities in the superconducting and normal 
states R s /R n s=s 1 above 100 cm -1 [16], in frame of the ISB model it means, that the high energy bound 
of the electron-boson interaction function Vt max < 100 cm -1 . The low energy part of the normal state 
reflectivity R(w), measured at T=40 K, was fitted (with the adopted fixed dj p i— 10000 cm -1 ) by (33) 
(see Fig. 3b), m* pi (w,T)/mi, is shown in the inset. One can see, that the optical mass at small u is 
huge itself and gives A > 100. This values points to heavy fermion or polaron models, rather than to 
the generally accepted weak or intermediate coupling scenario. On the other hand, the analysis of the 
normal state specific heat data gives an extremely small value of the coupling constant. As a result 
we concluded that Rb 3 Ceo is not a ISB metal and its properties should be treated in terms of a more 
complicated model which takes into account the self-energy and conductivity vertex corrections, the 
small value of the bandwidth AE «100 meV, possible influence of the electron-electron correlation, the 
strong anisotropy of the coupling function and the Fermi velocity, and etc. 

And the last but not least, if the data can be fitted by (33), that is by the formula with four parameters, 
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it is impossible to recover from this four parameters the detailed shape of the electron-boson interaction 
function a 2 F(cu). 

In conclusion, it is naturally to define the amount of the available information by the number of 
poles and zeros of the function which fits well the experimental curve or the number of the determinate 
Matsubara values. 



D. Temperature dependence of the optical mass and relaxation rate and its two-band ersatz 
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FIG. 4. Panels (a) and (b) show the frequency dependence of the optical mass and the optical scattering rate 
for T=100K, 200 K, 300 K. Panels (c) and (d) present the corresponding temperature dependence of m* pt (T)/mb 
and 7o P t(T) for lo = and lo — oo. All curves (solid lines) was calculated using model 2 spectral function. For 
comparison the normalised optical scattering rates 7*(o»,T) = 7 op t(c 1 j, T) / m* pt (uj , T)/nib are plotted by dotted 
lines in panel (a). The temperature dependence of the quantity o»( m *_i)/2 (see text) is shown by solid line on 
inset (c) together with the phenomenological line u> = 5T (dotted line). The common for r y op t(iJ = 0,T) and 
j op t(u = oo,T) asymptote 7 = 2-k\T is shown by dotted line in panel (d). 

In frame of the ISB model we assume, that all quasiparticles have the same Fermi velocity, but in 
real metals there is some anisotropy of vp. If the impurity scattering rate jimp is small in comparison 
with the complete optical scattering rate 7 opt and the electron-boson coupling is isotropic, we can use 
the average value < v 2 F > and keep the ISB model. On the other hand, if the impurity scattering, given 
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by the mean free path I, dominates, the quasiparticlcs from different parts of the Fermi surface will 
have different scattering rates ji m p,i = vf,z/1- Let us accept for simplicity, that we can divide the Fermi 
surface between two parts with the approximately uniform properties. The optical conductivity of this 
two-band system is the sum of two Drude terms 

a(io) = — — — ; h — — — — . (36) 

7l — ILJ 72 — ILU 

One can see, that Eq. (36) and Eq. (33) are equivalent and four parameters w 2 ( 1 ,uj^ 1 2 , 71, 72 in (36) 
connected with co^, A, B, C in (33) as follows 

2 2 2 

Upl = Upl.l + w pi,2 

c = ^pi,i72 +^ t|2 7i 

W pi,l +UJ pl,2 



A = 



7172(^,1+^,2) 
w p;,i72+^ i2 7i 



S = 7l +72- + -C (37) 

It means, that even if there is no electron-boson interaction in the two-band system, the frequency 
dependence of its optical conductivity (36) will reproduce a ISB model one (8) with some frequency 
dependent "mass" and "scattering rate". The multiband ersatz of the electron-boson effects could 
be revealed by the investigation of the temperature dependence of m* pt (w, T)/rrib and j opt (uj,T). It 
originates mainly from the presence of the Bose factor in Eqs.(5) (= sampling in (6)) and to a less 
degree from the Fermi distributions in Eqs. (5-9). 

In subsection II B the low (T — > 0) temperature behaviour of m* pt (u>)/m,b and j opt (uj) was discussed 
in detail. We found, that at a given quasipartical frequency u> in accord with the frequency conservation 
principle only virtual bosons with f2 < to can contribute to the relaxation. 

At high temperature the bosons lose their individuality. At T — > 00 the expression (5) can be 
substantially simplified, since the Bose factor coth(z/2T) w 2T/z dominates in comparison with the 
Fermi one tanh[(cj — z)/2T}. Keeping the leading term 2T/z we arrive at 

j -, n , x 7*m P , rr f dya 2 F(y) -y lmp + 2n\T 

Imuj{uj) ee 7e//(w) w + itT j = . (38) 

In Matsubara technique the same result can be obtained even more easily. At T — > 00 all kernels of 
the spectral function Afc, defined by Eq. (7), vanish as l/(< ^l oson > +47r 2 fc 2 T 2 ) for all k, except k=0. 
Therefore, the Eq. (6) for the renormalized frequency becomes trivial 

u> n = w n + 7„„ p /2 + n\T. (39) 

If one substitutes the obtained cu(ui) in Eq. (8) or u)(icj n ) into Eq. (9) after some algebra one arrives at 
the Drude type expression for the optical conductivity 

a(u) = ^ . (40) 

Note, that the relation -f opt = 2j e ff = — 2/mS is correct, if and only if the self energy S(w) = lo — Cj(lo) 
does not depend on frequency. Otherwise, the relation 

lopt{oj) = -2/m£(u;) (41) 

does not satisfy the Kramers-Kronig (K.-K.) relations, that is, it violates the causality. The proof is very 
simple. Due to the charge conservation law, when the photon has been absorbed, two quasiparticlcs, 
"the electron" and "the hole" come into being. The optical conductivity is a bosonic function defined 
by its values at v kl while the quasiparticle self energy T,(lo) is a fermionic function defined by the values 
at LJk- As a result, the right and left hand sides of Eq. (41) have different statistics, what is very 
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strange. In frame of the ISB model considered here the proof is a little bit more complicated. The 
optical conductivity itself and its inverse l/a(u>) are response functions. It means, that they have no 
poles on the upper half-plane of the complex frequency or, in other words, their real and imaginary parts 
are connected by Kramers-Kronig relations. The same is valid for the self energy. From Eq. (41) after 
K.-K. transformation one obtains 

m* opt (w) - 1 = 2(m* eff (u) - 1). (42) 

The optical and effective masses have the same asymptotic value 1+A at ui, T — > 0. Due to the presence 
of the factor of 2 in (42) it would mean, that 1=2. 
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FIG. 5. The frequency and temperature dependencies of the two-band optical mass, scattering rate, reflectivity 
and DC resistivity are shown by solid lines in panels a, b, c, d respectively. The corresponding ISB fitted 
counterparts are plotted by the dotted lines. 

Let us turn to the temperature dependence of the optical properties of ISB metals, which are shown 
in Fig. 4. All curves were calculated using model 2 spectral function and 7i mp = 0. Panels 4a presents 
the frequency dependence of the optical scattering rate 7 opi (u>,T) (solid lines) for T=100 K, 200 K, 
300K. The corresponding optical masses are shown in panel 4b. One can see, that the absolute value of 
-f op t{u,T) increases with T, but the frequency dependent component decreases. It leads to diminishing 
of m* opt {ui , T) / 'm& when T is rising. The temperature dependence of the minimum (lu = 0, solid line) 
and maximum (uj = oo, solid line) values of j opt (uj) are shown in panel 4d in comparison with their 
common asymptote 7 = 2-7rAT (dotted line). The corresponding function m* pt (w = 0,T)/rrib, named 
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sometimes 1 + X op t{T) is presented in panel 5c. The position of maximum of the small hump at low 
T is approximately equal to f2 m in/5, where fi m i n is the frequency position of the lowest boson peak in 
a 2 F(uj). At the same temperature the DC resistivity p(T) = 4irj op t(uj = 0,T)/u pl2 starts the linear 
rise. In subsection II B the simple rule for determination of the upper frequency bound of the spectral 
function Q max was proposed. Atwsi $l m ax the optical mass without unity (m* pt (Cl max ,T)/mb — 1) is 
approximately the half of its value at u = 0. The temperature dependence of the quantity W( m »_i)/2 
defined by 

2K pt (w = w (m ,_ 1)/2 ,T)/m 6 - 1)] = m* opt (uj = 0, T)/m b - 1 (43) 

is shown by solid line on inset in panel 4c in comparison with the phenomenological line u> = 5T. One 
can see, that 

W( m ._i)/ 2 « max(5T, Q max ) (44) 

The normalised optical scattering 7*(w) = 7 op t(w) / 'm* opt {uj ,T) / 'mb is shown in Fig(4a) (dotted lines). 
In our opinion [7] this quantity is far useless in comparison with j opt (uj) and m* opt {w, T)/rrib- With the 
exception of the low frequency region where its value is constant, 7*(w) is a quasilincar function ocwup 
to max(4Q maxl 2ir\T). Even the simple analysis by the favourite polynom Eq. (33) allows to determine 
the four parameters, while 7*(w) depends on two. It is interesting, that the Gilbert transformation from 
7*(w) is not equal to 0. In this connection we note that the use of the BSC normal state conductivity 
expression 

*2 

a{uj) = -, (45) 

v ; 4tt(7* -iu) 

is correct, only if 7* = const. Here uj*f = u pl / (1 + A) is the frequency-independent renormalised plasma 
frequency. 

Note, that according to the Allen formula (31), the optical relaxation is approximately linear with 
respect to the spectral function a 2 F(u>). It means, that the conversion from the low temperature scenario 
to the asymptotic high temperature one takes place for each boson with frequency fli separately at 
T w fi»/5, that is, when becomes of the order of the first Matsubara energy v\ = 2nT. 

Since the isotropic single wide band model is the subject of the present lecture, discussing the prop- 
erties of two-band systems we have to pay special attention to possible artefact results. Let us start our 
consideration from the temperature independent two-Drude terms simulator (36). Since it is equivalent 
to the favourite polynom (33), we can "extract the information about spectral function" including the 
temperature dependence of the coupling constant A(T). It is known, that if the material parameters are 
temperature independent, the optical scattering rate 7 op t(w,T) is a monotonically increasing function. 
Since the simulator's 7 op t(w) does not changes, when T is rising, one may conclude, that the coupling 
"constant" A(T) is a monotonically decreasing function with the high temperature asymptote 1/T. On 
the other hand, the analysis of the optical mass, which is usually a monotonically decreasing function, 
gives the opposite result: A(T) dramatically increases when temperature is growing. The study of the 
reflectivity R(u>,T) gives the third fitted function A(T). The reasonable way to resolve this conflict is 
the replacement of the model. If one restricts oneself to the consideration of the scattering rate only, 
the contradiction can be overlooked. The fact is that the strong temperature dependence of material 
parameters exists, for example, in heavy fermion systems. Another interesting scenario was realised in 
fullerides, where the DC resistivity (measured at the constant pressure) changes with temperature as T 2 
and the same DC resistivity (measured at the constant volume) follows conventional linear law p cx T, 
as it should be, since the theory assumes the V=const condition. 

The second example is the clean limit of anisotropic systems. As usually, one can divide the Fermi 
surface between two (for simplicity) parts with the approximately uniform properties. This disjoined 
representation of the Fermi surface is known as the multiband model. We accept also, that the intra- and 
intcrband interactions are governed by the electron-boson spectral functions having the same "model 
2" shape, but different coupling constants Ajj. The scales of transport and optical properties of the 
bands are fixed by the plasma frequencies Huj pl i = 47re 2 w|, i A r i (0)/3. In the normal state is possible to 
diagonalise the equations by the substitution Ai = A11 + A12 and A 2 = A21 +A 2 2- The resulting two-band 
optical conductivity is the sum of two ISB terms (8). For the model calculations the arbitrary values 
Ai = 2, A 2 = 1 and uo 2 t 1 /w 2 ; 2 =2 was chosen. 
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The frequency and temperature dependencies of the two-band optical mass, scattering rate, reflectivity 
and DC resistivity are shown by solid lines in Fig. (5)a, b, c, d respectively. The by-eye-analysis can 
not reveal the difference between clean limit two-band optical properties and the standard ISB ones. 
Fortunately the ISB model can pose several possibilities for a quantitative study. I simply fit the two- 
band curves, shown on panels (5)a, b, c, by the appropriate ISB ones (dotted lines), calculated using the 
same "model 2" spectral function. The coupling constant A was used as fitting parameter. One can see, 
that the temperature dependencies of m* pt (oj,T)/mb, 7opt(w,T) and R(uj 7 T) of the anisotropic system 
are yet weaker, that it should be according to the ISB rules. The fitted functions A(T) grow or drop, 
when T is rising, depending on which quantity has been chosen for the fit. 

For the few metals, where the mutual anisotropy of A(k) and of vp(k) was investigated, the following 
regularity was found. The quasiparticles with smaller vpfe) have larger A(k). Usually this information 
is unavailable. One can estimate the anisotropy, if one compares the material parameters obtained from 
the optical spectra with the ones revealed by the analysis of thermodynamic, or point-contact tunnelling 
data. The scale of the optical conductivity is fixed by uo pl oc v F N(0). The normal state specific heat 
data are proportional to (1 + \)N(0). The point-contact tunnelling spectra are weighted by the factor 
vp. As a result, if one treats optical spectra of anisotropic metals in terms of the ISB the model, one 
reveals the electron-boson coupling function a^ pt F(uj), which distinguish itself by smaller value of the 
coupling function A op t in comparison with lambdas, obtained by other physical methods. This difference 
is usually assigned to the kinetic coefficient (1-cos) [2] or, in another language, to the vertex corrections 
to the quasipartical Green function [1]. It is interesting, that if the system under consideration is 
completely isotropic, and a 2 F(w, k, k') do not depend on k and k', the vertex corrections vanish as soon 
as considered in this subsection usual anisotropy effects. What kind of corrections are more important 
in each particular case, first order vertex ones or discussed here zero order distortions caused by the 
interplay of w/(k) and A(k), it is the question for the gallop band structure calculations. 

In conclusion, the temperature dependence of the optical properties of ISB metals really exists and 
can be well quantitatively described. It originates mainly from the presence of the Bose factor in Eqs.(5) 
(= sampling in (6)) and to a lesser extent from the presence of the Fermi distributions in Eqs.(5-9). The 
principal deviation from the predicted frequency and temperature behaviour manifests that the chosen 
ISB model is not complete and has to be replaced. 



E. Optical mass and scattering rate in mid-infrared: artefacts generated by overlooking of 
interband transitions and underestimation of e^. 



The formal textbook relation between the conductivity <j(lu) and permeability e(uj) 

e( W ) = 1 + (46) 
to 

is not suitable for a practical use, since the entering Eq. (46) <t(oj) reflects the total response of all 
excitations. 

The simplest reasonable decomposition of (46), valid over the spectral region up to frequencies just 
above plasma edge, is 

i \ i \ i \ 4:Tria(Lu) 

^V^V = Winter {(d) ~t~ tphonon 

(u + — , 47 

where €i n ter and e p honon are the contributions of interband transition(s) and direct IR active phonons 
centred below plasma edge. The response of high energy excitations and wide interbands, overlapping 
the plasma edge is presented by the complex constant eoo. cr(w) is the quasiparticle conductivity, defined 
by Eqs.(8,10). 

The direct phonon band contribution will not be considered, since it is usually weak and plays some 
role only in the selected compounds containing light elements and having small value of uj p i, for example, 
in doped fullerencs. 

The used here model interband contribution 

einter M = fi2 ^jf"^ (48) 
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is the simple Lorentzian with u> p i^ n t e r =16000 cm" 1 , S!o=4000 cm -1 , and 7i„ ter =6000 cm' 1 . The over- 
looking of the interband contribution was simulated as follows. The sum of two terms given by Eqs. 
(48) and (8) was treated as the quasiparticlc response in terms of the optical mass and scattering rate 

V(^) - ium* op MT)/m b = — ^j^— (49) 

For simplicity I set lo p i = uj p isnter and following to 

The results are shown in Figs. 6a, 6b and 6c by solid lines 

j opt (uj,T) - iuim* opt {u),T)/m b = 

and the chi-by-eye [19] fit by the formula (33) are plotted in Figs. 6a, 6b and 6c by dotted and dashed 
lines respectively. 



. The "true" masses and scattering rates 
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FIG. 6. The artefact frequency dependencies of the optical masses and scattering rates (insets) generated by 
the overlooking of interband transition (panels a, b, c) and the underestimation of e x (panel d) are shown by 
solid lines. The true, unpertubed frequency dependencies of m* opt (uj)/m b and 70^4(1^1) are plotted by dotted lines. 
The dashed lines in panels a, b, c are the chi-by-eye fitted curves ( by Eq. 32). For parameters of the model 
interband transition and trial ISB quasiparticle optical conductivities see text. 



Note, that due to the frequency averaging, the quasipartical optical conductivity (8,10) is a slowly 
varying function in the mid-infrared region. In contrast, the interband permeability near transition edge 
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(here, near the center of the Lorentzian) can changes rapidly. The second important feature of ej n t er (o;) 
is its very weak temperature dependence. 

The mixture of the standard Drude term (7=1000 cm -1 ) with the Lorentzian (48) is shown in panel 
6a. This simple artefact can be easily recovered, since it is temperature independent and the optical 
mass changes too quickly. 

The second and third examples are the half- by- half blends of the "model 2" ISB conductivities (T=100 
K, panel 6b;T=300 K, panel 6c) with the same ei n ter{uj) Eq. (48). It is more complicated case, since the 
nonmonotonic behaviour could be screened by forthcoming transitions. In this case detailed quantitative 
analysis of both the frequency and the temperature dependencies is needed. One has to pay the special 
attention to the temperature dependencies of m* opt {ui — > 0,T)/m,b, U( m *-i)/2 an d 7 op t(o; ~ 2ir\T,T). 

If the interband transitions are wide spectral bands without unqualified features, their contribution 
to the permeability is usually included in . The effects caused by the underestimation of its real part 
was considered by us in [7]. Here we study the problems produced by the use of Eq. 46 and by the 
underestimation of the imaginary part of . Since the last quantity can be approximated by the simple 
constant only in the limited frequency range, for the sake of the analytical correct the model function 
was chosen in the form 



7771(^00(0;)) 



, /w- 2000\ , / 50000 
tanh — — — — + tanh — t^^t — 
V 5000 / V 30000 



at oj > 0, 



7m(<5e 00 (-w)) = -~Im{8e oa {uj)) 1 
7777(^00(0)) =0, 

Re(Se^)) = 2.5 - ^ f ^(S^x)) 

In the mid infrared [dc^u))) w 4 + 2i. 

The model frequency dependencies of the optical masses and scattering rates 

lopt (u,T)-icjm* opt (u,T)/m b = - p [ (52) 

47T(7(CJ) — lUJO^oo) 

are shown in panel 6d for 5{eoc) = S(e oc (uj)) — 1 (solid lines), iffjo) = S(e 00 (cj)) — 4 (dashed-dotted 
lines), and <5(eoo) = 2? (dashed lines). The unperturbed 7 OJ , t (w,T) and m* pt (u,T)/mb ("model 2" 
spectral function, A=2, T=300 K) are plotted by dotted lines for comparison. As usually, the masses do 
not lie, and manifest that the "unusual" frequency dependence of the scattering rate is no more than 
the artefact of the chosen interpretation procedure. 

The early band structure calculations [27] predicted, that the interband contribution into the in-plane 
HTSC permeability is a real constant between 10 and 15 at u — > 0, and becomes complex for u > OAeV. 
In the vicinity of the plasma edge its value is 4+2L This information was used by us for the successful 
description of the normal state optical properties of HTSC materials both in FIR and MIR spectral 
regions [7,15]. 

In conclusion, the overlooked interband transitions produces the "unusual" frequency dependence of 
the optical relaxation, which can be nevertheless analysed properly within the ISB model. 



F. Why is the first derivative from the experimental curve better than the second one? 

Let us rewrite the valid at T=0 Eqs. (29,30) in the differential form 

JW-j.tZtM (53) 

At a glance both Eqs. (53) and (54) look nice and simple. But if one attempts to apply them to an 
experimental data, the distinction becomes evident. The point is that the experimental data are known 
only approximately. A finite value of the signal to noise ratio (S/N) is the reason why the Eq. (54) 
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can not be applied to the FIR data. Typical value of the reflectivity S/N do not exceed 50, the multi- 
reflection techniques [3] allow to reach S/N=2000. Let us adopt for simplicity, that f(u>) = u"f opt {uo) is 
sampled in equidistant points u>i, i = 1, M and all values /(w,) have the same errors 



/(Wi) (S 



(55) 



Note, that in the normal state at T = (when the application of (53) and (54) is correct) both 
f(uj) = ujj opt (uj) and its first derivative f'(co) are monotonically increasing functions, if a metal under 
consideration is the ISB one. 

The difference f(cj i+ i) — has the order of f(ui)/M. Its error is the ^/(5/ 2 (o> i+ i) + Sf 2 (u)i) ~ 

y/25f(ui). As a result the first derivative f'(uj) has the signal to noise ratio 



(s\ _ rim) _ /(^ +1 )-/(u> t ) ^ i rs\ 



(56) 



reduced by the factor K\ = \/2M in comparison with the original (S/N) in f(w). Similar the second 
derivative /"(w) has the signal to noise ratio 



s_ \ f"( Ui ) _ (,r^ +1 )-,f(^))' i (s_ 

NJ f » Sf»(wi) y/SfVfa+J+SfVfa) ~2M*\Nj f 



(57) 



reduced by the factor K 2 = 2M 2 . 

In subsection II C the four parameters formula (33) was promoted as the universal fitting tool. M = 4 
means K 2 = 32. In another words, the normal state 7 opt (w) has to have the original signal to noise ratio 
(S/N)=100, in order to define four points in its second derivative 7o P ^(w) with accuracy of 30%. 

As for Ki, this quantity grows reasonably. Of course, one can not measure the effective scattering 
rate 7 e //(<*>) and use the formula (53), nevertheless the title question is worth-while. The point is that 
the first derivatives from the superconducting state optical scattering rate j' s pt (uj) and the reflectivity 
R'{io) reproduce the input spectral function a 2 F{uo) [2,3,20,21]. The well structured spectral function 
of lead obtained from the analysis of the first derivative of A s (w) — A n (u) looks convincing (at least for 
me) due to the reported value of S/N=6000. 

Analysing Eq. (54) we arrive at the same conclusions as discussed above. The normal state optical 
conductivity provides us with restricted information similarly as other physical methods. The only way 
to extract the all required material parameters values is the analysis of different measurements in frame 
of the same model. 



III. SUPERCONDUCTING STATE OPTICAL CONDUCTIVITY 



A. The inversion procedure 



1. Adaptivity, what does it mean? 

Let us consider from a general point of view the problems, which one has to take in mind, when one 
try to reconstruct a spectral function a 2 F{uj) from experimental reflectivity data. This analysis should 
be performed in frame of a model which is based on an appropriate formalism. Thus, it means, that the 
results obtained arc model dependent. The ideal, but rare situation occur, when the exact analytical 
solution of the inverse problem is available 

a 2 F(tj) = f(R(uj),uj) (58) 

where / is a known function, which could depend on a 2 F(uu). In the last case the numerical solution 
can be found by iterations. 

But what should we do, if such a function / is not available by reason of the complexity of the 
direct formalism. In this case we shell use the numerical adaptive strategy, based on the substitution 
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of equations of the direct formalism by the adaptive formula(s), suitable for the numerical solution of 
ill-posed problems. Since it is a strategy, let us consider how does the adaptivity work. 

Generally speaking, one can use any adaptive formula. For example, one can consider the a 2 F(oj) as 
an input function and the reflectivity as an output function for some formalism. Then one could write 

R(lu) = J dzK(uj- z)a 2 F(z), (59) 

or in the Fourier space 

T{R{x)) = F(K{x))F(a 2 F(x)), (60) 

where T(f(x)) is the Fourier image of the function / and K(x) is an for the time being unknown kernel. 
K(x) could be determined as follows. 

• At first, one has to substitute a trial spectral function ct 2 rial F{u) into the equations of the main 
formalism (Eliashberg equations + Nam-Rainer formula) and obtains a trial reflectivity R tr iai(u). 

• Next, one should substitute a 2 rial F(uj) and the calculated Rtriai{u) into the adaptive formula (60) 
and solve Eq. (60) for Ki(w). The subscript 1 means "the first iteration". 



Finally, one has to substitute the now known kernel Ki(w) and the experimental R exp (uj) into the 
same adaptive formula (60), and solve it for a\F(w). 



So, on each iteration we use the model formalism ones and the adaptive formula twice. At first, one 
have to calculate Rk-i{uu) using a\_ 1 F(uo), as input function. At second, one should find the adaptive 
function(s) 

T{Kk{x)) = H*UF(x)y (61) 



and finally calculate the spectral function 



F{R exp (x)) 



" "TOT (62) 

The presented "trick" is quite general. One could replace analogously, for example, R(u>) in (59-62) 
by the tunnelling density of states N(u>) or the ARPES spectral function A(u), and obtain adaptive 
methods for the tunnelling and electronic spectroscopies respectively, which is expected can work, with 
any formalism. 

As mentioned above, the experimental data are known only approximately, therefore the inverse 
problem is ill-posed. Up to now the solution of ill-conditioned problems is the art more than the routine 
procedure. Since one can use any adaptive formula(s), it will be pragmatic for him/her/AI to choose it 
from the short (alas!) list of well investigated by numerical mathematics equations (sec, for example, 
sec. 18.4-18.7 in [19] and references within). 

On this understanding the discussed above convolution adaptive formula is one of the best possible 
selections, since the Fourier transformation is linear and unitary. 



2. The inversion method for the s-wave ISB reflectivity. 



The inverse problem for the superconducting state optical absorptivity of Pb was successfully solved 
by B. Farnworth and T. Timusk [3] within the weak coupling ISB model, based on the theory proposed 
by P.B. Allen [2]. The difference between weak coupling theoretical calculations and experimental data 
[20] was analysed in detail in [2]. It was pointed, that: 

• The first derivative from the superconducting state optical scattering rate contains the component, 
proportional to the input spectral function a 2 F(uj) due to the influence of coherent factors (11). 

• The weight w/2Ao of this component is strongly underestimated. 
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• If one substitutes the solutions of (2) into the dominators of (10) and keeps BSC values cD = ui and 
A = A for the coherent factors (11), one obtains the optical relaxation, close to the one, given 
by the weak coupling Allen theory. 

Below we will see that ISB model based on the Nam formalism have qualitatively the same properties. 
There are several reasons, why the following convolution adaptive formula 

2^1e{uj) = 2ir J dzK{u - z)a 2 F(z)6(z) (63) 
9(x) = 1 if x > 0, else 6{x) = 

is suitable. At first, we segregate the ill-posed part of the inverse s-wave problem (taking of the first 
derivative) into the separate step. At second, the low and high frequency values of 7 p t (w) coincide (=0), 
it is important for the FFT algorithm. At third, the obtained kernel K (u>) reminds us the instrumental 
resolution function of a badly adjusted spectrometer, allowing to profit by the applied spectroscopy 
experience. If one thinks in terms of the transfer functions, Laplace transformation adaptive formula 
could be more suitable. 

The equation (63) was solved with respect to K (ui) using optical scattering rates calculated in strong 
(solid line, Fig. 7b) and weak (dashed line, Fig. 7b) coupling approximations and input "model 2" 
a 2 F(uj). The results are plotted by the same type lines in Fig. 7a. No doubt, all features of the kernel 
are important for the recovery of the precise peak amplitudes and energy positions. Nevertheless, the 
general shape and main physics of j^ pt (uj) are defined by the master peak, located just above ui — 2A . 
The arrow, pointing in Fig. 7a the frequency position of 2Ao, is located within this peak. The master 
peak income is illustrated by dotted line in Fig. 7b. It was calculated as follows. 
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FIG. 7. Panel (a) shows the frequency dependence of the superconducting state kernels K(co) (solid and dashed 
lines) in comparison with the normal state one (dotted line) . The source superconducting state optical scattering 
rates, calculated in the strong (solid line), and weak (dashed line) coupling approximations are shown in panel 
(b). The convolution of the input spectral function with the master peak (=" coherent" relaxation) is drown in 
panel (b) by the dotted line. 



• All collateral maxima and minima of the shown by solid line K{oj) were rejected and artificial kernel 
K co h(uj) was set to zero anywhere besides the master peak frequency region, where it coincides 
with K(w). 

• K co h(uj) was convoluted with the input spectral function a 2 F(uj) and integrated over co 

• Since the discussed phenomenon is governed by the coherence factors, the resulting "scattering 
rate" (dotted line, Fig. 7b) is labelled as coherent relaxation. 
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• The simple frequency dependence of the coherent relaxation, reminding us the properties of the 
quasiparticle effective scattering rate j e ff(u>), allow to enunciate the Visual Accessibility procedure, 
described in forthcoming section. 

The normal state kernel (multiplied on 3) is shown by the dotted line on Fig. 7a for comparison. The 
pessimistic conclusion, made in the previous chapter, can be couched in current terms: the (induced by 
the frequency and temperature averaging) width of K(ui) is too large to resolve the fine structure of the 
spectral function. 

The domination of the master peak makes the choice of the trial a 2 F(u>) on the first iteration easy. It 
can be the normalised 7^(0;) or R' s (ui). Note, that the master peak is asymmetric. The final frequency 
resolution is defined by the small value of its left half-width. 

The self-consistent determination of the value of the superconducting gap is theoretically possible in 
frame of adaptive approach, but strongly not recommended. Following this way one transforms the 
weakly nonlinear task to the strongly nonlinear problem. For the case of ill-conditioned equations it is 
the wrong way. The value of Ao has to be known a priori. The related issue is the difference between 
a 2 F(uj) and a 2 r F(u). Substituted to the Eliashberg equations transport spectral function a 2 r F(uj) lands 
us in difficulty, since the obtained value of the gap can be lower than the known A . In my opinion, for 
the sake of simplicity in this case one might sacrifice the unknown (from the first principle calculations) 
Coulomb pseudopotcntial and use fi* as a simple fitting parameter. 

In conclusion, the adaptive method of the exact numerical inversion of the s-wave superconducting 
state Nam equations is presented. In my opinion, one day the FIR spectroscopy will provide not an 
alternative to the tunnelling, but a powerful independent tool for the determination of a 2 F(uj) in s-wave 
superconductors. 



B. Rules of Visual Accessibility procedure 



The frequency dependence of the superconducting state ISB reflectivity R s (co) in the local limit has 
the distinctive forms. This features are determined by the physical properties of the optical scattering 
rate j^ pt (uj). In this section we consider our simplified approach (Visual Accessibility) which allow us 
to analyse "by eye" the experimental reflectivity curve and get qualitatively the same results as the 
sophisticated method, described in the previous section, able to obtain quantitatively. 



1. Optical scattering rate and reflectivity 

As far as the physics is determined by the "f op t((jj), but the measured quantity is R(w), let us derive 
a simple relation between them. In view of a high value of the permeability e(u>) 3> 1 one can expand 
the formula 



R(oj) = 



1 + VeR 



(64) 



in powers of l/^/e(uj). In view of Eq. (22), substituting 1 + A instead of m* opt /m,b, one arrives the 
following approximate relation 



R(lo) w 1 - 



_27optH_ 
Wpi\/l + A 



(65) 



valid at u> < Vl max + 2A both in the normal (A=0) and superconducting (A=Ao) states, except for 
small frequencies w < 7 op t(0)/(l + A), where Hagen-Rubens approximation is suitable. 

Fig. 8a shows the frequency dependence of the superconducting state reflectivity i? s (w) (T=10 K, 
solid line) in comparison with the normal state reflectivity R n (u)) (T=100K, dashed line). Both curves 
have been calculated using the "model 2" spectral function. The frequency dependencies of the optical 
scattering rates jg pt (u>) and 7™ pi (w) are shown in the inset. One observes that R(u>) looks like the mirror 
image of the appropriate 7 op t(w), reproducing all important features. 
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FIG. 8. Panel (a) shows the frequency dependence of the superconducting state reflectivity R s (u>) (T=10 K, 
solid line) in comparison with the normal state reflectivity R n (co) (T=I00K, dashed line). Both curves were 
calculated using "model 2" spectral function. The frequency dependencies of the optical scattering rates 7„ P ((w) 
and 7o P t(w) arc shown in the inset. Panel (b) shows the frequency dependencies of the first derivatives from 
the optical scattering rate j' s pt (u))/2n and from the normalised reflectivity R'{u))u p i(l + Ao)°' 5 /47T in comparison 
with the shifted on 2Ao input spectral function a 2 F(cj — 2A). One can see, that due to domination of the 
coherent relaxation in the frequency region between 2Ao and 2Ao + n max , i?'(w)w p i(l + Ao) /4-k and 7op t (o))/27r 
reproduces well the peaks of the spectral function. The important features are indicated in panel (a). One can 
see, that the absorption onset at 2Ao, the frequency dependence of the coherent relaxation between 2Ao and 
2Ao + tlmax, and its upper energy bound at 2Ao + Qmax are visually accessible. The term "integrated spectral 
function" means, that the first derivative over u> from this part of the spectrum reproduces a 2 F(cj — 2Ao). 



2. The visual accessibility criteria for s-wave superconductors 

Let us summarise our knowledge about the s-wave superconducting state optical relaxation rate and 
the local (London) limit reflectivity (see Fig. 8). The following features are visually accessible. 

• The doubled gap value 2A coincides with the absorption edge. If the u < 2A , the reflectivity 
R(w)=l and 7$pt(w)=0. 

• The impurity scattering gives a sharp decrease of the reflectivity just above 2A . The magnitude 
of the "jump" is proportional to 7, mp . -fg pt (uj) jumps approaching 7 im? , from above. 

• In superconducting state the coherent relaxation dominates. Hence, the first derivatives from 
7op t (u;) and from R(w) reproduce the shape of the spectral function (see Fig. 8b). 

• Since the frequency dependence of the coherent relaxation rate grows faster, than the frequency 
averaged normal state one, and is fulfilled its maximum value at the frequency to — 2A.Q "f^max? 

at 

this point the ratio of the -f^ pt (uj) / -f^ pt (uj) has a maximum, while the reflectivity ratio R s {u)/R n {w) 
exhibits a minimum. 



C. Visual accessibility and separable d-wave model 



In this section we consider briefly, how the procedure described in section III B 2 should be corrected 
in order to make the main features of the rf-wave ISB reflectivity visually accessible. 
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1. Separable d-wave model 
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FIG. 9. Panels (a) shows the frequency dependency of the reflectivity R s (u>) (dashed lines) calculated within 
the separable d-wave model in comparison with the experimental data (solid line, Schutzmann et al Phys. Rev. 
B46 512 (1992)). The spectral functions aoF(u>) are shown in the insets. The coupling constants were chosen 
to obtain the value of the tunnelling gap A tun =200 cm -1 . The values of the plasma frequencies lo p i and the 
low-frequency values of the non-metallic part of dielectric permeability e x were obtained by the least-square 
fit. Panel (b) shows the frequency dependence of the first derivative from the calculated optical scattering rate 
7opt(w) (dotted line) in comparison with the shifted on 2A opt input spectral function a^F{ijj — 2A op t)(solid 
line). The convolution of the o^F{uj — 2A opt ) with the Gauss distribution function (shown in inset (b)) is 
plotted by the dashed line. Panel (c) shows the frequency dependence of the calculated tunnelling density of 
states N(w), the frequency positions of the tunnelling gap Atun and the optical gap A op t. Panel (d) shows 
the frequency dependencies of the reflectivity R s (ai) (dashed lines), calculated using caudate spectral function 
(inset), in comparison with the experimental data. 

Formally speaking, the frequency dependencies of the isotropic and the anisotropic parts of the spectral 
function are not obliged to coincide. At the same time, if the shapes differ substantially, the density of 
states (67) looks strange. So, for the sake of tunnelling, the model spectral function was chosen in the 
simple form 

a 2 F{uj, if) = alF(uj)(\ + 2gcos(2f)cos((p / )) (66) 

If one substitutes the function (66) into the s + id ISB Eliashberg equations, one finds that at g > 1 the 
d-wave component entirely dominates. Hence, instead of the sophisticated s + id approach, we can use 
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for our analysis the pure d-wave model described in detail in [31]. For the sake of simplicity we restrict 
ourselves here to the consideration of the clean limit 7i m?J = and set the Coulomb pseudopotential 
//*(</?, tp') to zero. The value of the parameter g was chosen little a bit larger than unity, since the normal 
state properties (defined in frame of this model by a^F{uj)) manifest the strong renormalization due to 
the phonon-electron interaction. 

It is interesting, that it is possible after minor rehash to use the available s-wave computational 
programs for the strong coupling d-wave calculations. The clean limit version of the separable d-wave 
model keeps the convolution structure of the set (2). One has to substitute the angle averaged density 
of states 

/ &(y) \ 1 f u>(y)dtp 



JO 



\yJu,(y)Z-A*(y)/ ^ * Jo yj&(y) - 2A 2 { y )cos 2 {2^) 
instead of N(ui) (3) in (2), and the angle averaged density of pairs 

/ My) \ 1 r 2A(y)cos 2 (2(p)dip 



^{y) 2 -A 2 {y)l ^ n Jo ^ 2 ( y )-2A 2 ( y )cos 2 (2<p) 



(68) 



instead of D(u) (4) in (2). Similarly, the conductivity can be calculated using the Nam formula (10). 
For this purpose one should substitute A<j(x) = \/2A(x) cos (2tp) instead of the A s (uj) into the element 
of integration in (10) and averages the obtained cr(u), <p) over the angle tp. 

Fig. 9a shows the model reflectivity R s (w), calculated in the frame of the separable d-wave model 
(dashed line) in comparison with the experimental data [5] (solid line). The spectral function ciqF{uj) 
is shown in the inset in Fig. 9(a), its coupling constant A was chosen to reproduce the value of the 
tunnelling gap Ai„„=200 cm -1 derived from the tunnelling spectra [32]. Note, that it was the same 
"model 2" [7] curve, used above, but without any high energy peak. The value of the plasma frequency 
ujpi and low frequency value of the non-metallic part of dielectric permeability e M were obtained by the 
least-square fit. 

Fig. 9b shows the first derivative from the d-wave optical scattering rate J opt (u>) (dotted line) in 
comparison with the shifted on 2A opt input spectral function ocqF(uj — 2A opt )(solid line). Here we define 
the optical gap A opt « A tun /V2 as the physical quantity, with defines the shift of the spectral function 
alF{uj) in FIR and tunnelling spectroscopy. The distinction between A opt and A 

tun &riS6S due to the 

angle averaging (< ... > v ) both in the Eliashberg equations and in the conductivity calculation. The 
tunnelling gap A tun (see Fig. 9c) is the maximum value of the d-w&ve gap 

A d ((p) = V2Acos(2ip) (69) 

and the optical one is its mean square angle average. Note, that the dip in the density of states N(w) (see 
Fig. 9c) arises at A opt + £l ma x and the electron-boson region starts in N(u>) at A opt < A tun . It means, 
that the low-energy excitations (if they contribute to the spectral function) will degrade the sharp peak 
in the tunnelling density of states. 

If one substitutes J opt (u>) in Eq. 63, one finds that coherent relaxation still dominates, but the main 
peak is substantially wider and can be approximated by a Gaussian function Qexp{— uj 2 / A 2 opt ) ((Fig. 
9b, inset). The function 

C(w) - a 2 Q F(iu - 2A opt ) ® \ exp (-£-) (70) 

y/TT A opt \ A opt J 

is plotted by the dashed line in the Fig. 9b and agrees well with the 7 opt (u>) / '2ir . Here ® means a 
convolution. Note, that the coincidence of the prefactor Q with the l/y^Ao^ is come-by-chance. 

In my opinion, the additional ill-posed step, the real deconvolution of the Gaussian can not be per- 
formed using the data with reasonable value of (S/N), since the spectral function in these compounds 
have half- width of the order of A opt . On the other side, if one (by chance) will find that 7 opt (uj) or R'(u>) 
are structured, it means, that the ISB separable d-wave model is not suitable for HTSC and have to be 
replaced. 
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Strictly speaking, the d-wavc reflectivity is less than unity at any frequency. At the same time, there 
is an onset of the absorption at some wa, where R(o>) becomes visually different from unity. Since the 
convolution C(to) is wider than the input spectral function a 2 F(uj), the position of the onset u>a < 2A opt . 
For the model spectral function shown in the inset in Fig. 9a u>a ~ l-4A opt A tlm . Generally speaking, 
the selected boson modes could be responsible for the d-wave pairing. In this case the position of the 
absorption onset will be between f2 m j n + A opt and £l m i n + A tun , where TO j„ is the low frequency bound of 
the electron-boson interaction function. Similarly, the point where the coherent relaxation rate reaches 
its maximum, will be shifted by A opt towards the high frequencies in comparison with its position in 
the s-wave superconductors. As a result, the reflectivity ratio R s (uj) / R n (ui) exhibits a minimum and 
the scattering rate ratio j^ pt (uj) / j" pt (uj) has a maximum approximately at fl m ax + 2A tun . 

2. The visual accessibility criteria for d-wave superconductors 

In ISB s-wave metals at the same frequency u> — Aq the density of states exhibits three important 
features. At first, in this point N(u>) has a maximum. At second, since the averaged over k value of 
Ao coincides in this case with itself, the density of states at lu — Aq changes most rapidly. At third, 
for lu < Ao, the value of N(ui) is equal to zero. In rf-wave superconductors this structural feature has a 
different energy position. One could denote 

• the position of the maximum in the density of states Atun as the tunnelling gap; 

• the value of the shift 2A ov t of the image of the spectral function in the first derivative from the 
optical relaxation rate j opt (co) as the optical gap , 

• the frequency to a where R s (w) becomes visually distinct from unity as the onset of the absorption 
or as the absorption gap. 

The approximate relations between the gaps are A tun ~ ~ V^A opt . Similar to the description of 
the visually accessible features made in subsection III B 2 for a s-wave metals, we can write the same 
criteria for the separable d-wave model. Note, that in comparison with the s-wave case, all presented 
below itemised statements are valid approximately only. 

• The single gap value A tun coincides with the absorption edge. If the u) < A tun , the reflectivity 
R(w) wl and Yo P M «0. 

• The impurity scattering does not give the sharp decreasing of the reflectivity just above the A tun . 
Instead, it decreases the gap itself. 

• The coherent electron-boson scattering dominates in the superconducting state. The first deriva- 
tives from the Y opt (co) and from the R(w) reproduce the spectral function shifted on 2A opt and 
convoluted with the Gauss distribution function exp{— uj 2 / A 2 opt ) (see Fig. 9b). 

• Since the frequency dependence of the coherent relaxation rate grows faster than the normal state 
one, and reaches its maximum value when the frequency u w 2A ttm + Ct max , at this point the ratio 
of the ll p t{ w ) Hopt{ w ) wu l exhibit a maximum and reflectivities ratio R s {u)/R n {w) will exhibit a 
minimum. 

Finally, let us make the simple acceptance test for the presented above visual accessibility procedure. 
The test spectral function has a long history. At that old time, when people did not distinguish between 
the effective scattering rate j e ff(w) and the optical one 7 pt(w), Collins et al [22] interpreted the high 
energy asymptotic behaviour of the 7 op t(u;) (see Fig. 2) in terms of its step-like effective counterpart 
j e ff(uj) (5, 30). As the result, they arrived at the spectral function having a long tail up to the very 
high frequency (see Fig. 9d, inset). Nowadays, similar spectral functions occasionally arise in different 
models, since this shape suggests the spin fluctuations spectrum reported by neutron spectroscopy. 

The calculations was performed in the same fashion as for model spectral function shown in Fig. 9a. 
The coupling constant was chosen to reproduce the value of the tunnelling gap A tM „=200 cm -1 . The 
values of the plasma frequency lu p i and low-frequency value of the non-metallic part of the dielectric 
permeability e M were obtained by the least-square fit. The peak position had the value 330 cm -1 which 
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was close to the one in the original paper [22]. I choose the power law 1/w 2 5 for the tail and the high 
energy cut-off 1600 cm -1 . 

Let us compare the calculated curve and the experimental one. Since the low frequency bound in the 
model spectrum f2 m j„ w 200 cm -1 , the absorption edge takes place at Sl m in + A tu „ in accord with the 
results discussed in [26]. At second, in the region between 400 cm -1 and 800 cm -1 , where the reflectivity 
"integrates" the spectral function, there is a systematic shift on 140 cm -1 . For the compensation of this 
shift one has to decrease the value of tunnelling gap by a. factor of two The most important discrepancy 
is connected with the upper frequency bound of the single-particle relaxation, that is, with the Q ma x- 
The reflectivity ratio R s /R n has a minimum at 1400 cm -1 in comparison with the 900 cm -1 in the 
experimental one. 

In conclusion, the separable rf-wave model reproduce all important features of the traditional ISB 
s-wave one with the single exception: the resulting images of the spectral function are the convolution 
of the input electron-boson interaction function and the Gauss distribution exp(— uj 2 /A 2 ). 
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